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Strongly coupled field theories with gravity duals can be placed at finite density in two ways: 
electric field flux emanating from behind a horizon, or bulk charged fields outside of the horizon 
that explicitly source the density. We discuss field-theoretical observables that are sensitive to this 
distinction. If the charged fields are fermionic, we discuss a modified Luttinger's theorem that 
holds for holographic systems, in which the sum of boundary theory Fermi surfaces counts only 
the charge outside of the horizon. If the charged fields are bosonic, we show that the the resulting 
superfluid phase may be characterized by the coefficient of the transverse Magnus force on a moving 
superfluid vortex, which again is sensitive only to the charge outside of the horizon. For holographic 
systems these observables provide a field-theoretical way to distinguish how much charge is held 
by a dual horizon, but they may be useful in more general contexts as measures of deconfined (i.e. 
"fractionalized" ) charge degrees of freedom. 



CONTENTS 

1. Introduction 

A. Luttinger deficit 

B. Deficit in the Magnus force 

II. Setup 

III. Bulk fermions and Luttinger's theorem 

A. Proof 

B. Applications 

1. Confining geometries 

2. Electron stars 

3. Reissner-Nordstrom black hole 

4. Dilatonic systems 

IV. Vortices in holographic superfluids 

A. Background on holographic vortices 

B. Berry phase 

C. Magnus force 

D. Applications 

1. Finite temperature superfluids 

2. Confining superfluids 

3. Speculations at zero temperature 

V. Conclusion 



4 
4 
7 
7 
7 



9 

10 
11 
11 
13 
13 
14 
14 

15 

16 



A. Technical details of spinors 

1. Green's function evaluation of charge density 16 

2. Finite temperature Luttinger theorem 17 

B. Technical details of vortices 19 
References 20 

I. INTRODUCTION 

This paper will be concerned with holographic descrip- 
tions of compressible phases of matter. Recall that a 
compressible phase is conveniently defined to be one in 



which a [/(I) charge density p is a continuous function 
of various parameters such as the chemical potential /i 
and the temperature T. When the U(l) symmetry is 
unbroken, as emphasized recently in [1], all known field- 
theoretical examples invariably involve Fermi surfaces 
that carry the U{1) charge, and one expects Luttinger's 
theorem, which relates the charge density to the volume 
enclosed by Fermi surfaces, to apply. 

On the other hand, if the U{1) symmetry is broken, 
the system is in a superfluid state. While Luttinger's 
theorem no longer applies, there nevertheless does exist 
a probe of the charge density, which is provided by a vor- 
tex. If one moves a vortex in a closed loop the many-body 
wavefunction acquires a Berry phase; relatedly, a moving 
vortex in a superfluid feels a transverse force called the 
Magnus force. As we explain later, this Berry phase and 
Magnus force can be viewed as probes of the charge den- 
sity. 

The goal of this paper is to examine the status of these 
relations in holographic systems. We find that in both 
situations (whether the global U{1) is broken or not), 
these relations develop anomalous contributions which 
are associated with horizon degrees of freedom in the 
bulk, and which can in turn be interpreted in the bound- 
ary theory as deconfined (or "fractionalized") degrees of 
freedom with respect to the gauge group of the boundary 
theory. 



A. Luttinger deficit 

A fundamental result of Fermi liquid theory is Lut- 
tinger's theorem [2], which in its original form states 
that to all orders in perturbation theory the volume V 
enclosed in a Fermi surface of electron is equal to the 
charge density p: 



-V = p. 



(1.1) 



The original proof of Luttinger and Ward was based 
on perturbation theory. More recently an elegant non- 
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perturbative proof was given by Oshikawa in [3], which 
highhghted the topological nature of the theorem. Os- 
hikawa's proof also admits a generalization to exotic 
phases of matter in which the electron fractionalizes into 
constituent particles with different quantum numbers. 
Such a fractionalization is accompanied by an emergent 
gauge symmetry. It was found that in such a case Lut- 
tinger's theorem can be violated, with the deficit being 
related to the momentum carried by a topological exci- 
tation of the emergent gauge theory [4-6]. 

We now turn to holography. The field theories which 
are known to admit gravity duals all involve non-Abelian 
gauge theories with various matter degrees of freedom [7] . 
In such a theory we will define the location of a Fermi 
surface as a shell in momentum space with a radius kp 
which satisfies 

G-\LJ = 0,kF) = (1.2) 

where G is the retarded Green function of any gauge in- 
variant fermionic operator. The generalization of the 
Luttinger theorem (1.1) to such theories would then 
naively amount to to 

(27r)^-i J2<li^i = P (1-3) 

where the sum is over all Fermi surfaces of charged 
fermionic operators of charge qi. 

Now recall that the charge density p of the boundary 
theory is mapped to the boundary value of the radial 
electric field in the gravitational bulk, i.e. the gauge- 
gravitational dictionary gives us 

p=^V^F'^{r^oc) . (1.4) 

Here is a bulk gauge coupling and F^^ is the radial 
electric field component of the bulk U{1) gauge field cor- 
responding the boundary current. Thus if we are study- 
ing a field theory state with a nonzero p, something must 
create a nonzero electric field in the bulk. There are two 
known ways to implement this: 

1. One may study a system with a charged black 
hole horizon in the bulk. In this case the source 
of the electric field is actually hidden behind the 
horizon. Nowhere explicitly in the bulk geometry 
does one actually see any charged matter, and the 
electric field flux is thus constant throughout the 
geometry. One example of this is the Reissner- 
Nordstrom black hole [8] , which is characterized by 
its AdS2 X infrared geometry at zero tem- 
peratures. A different example is provided by the 
dilatonic systems studied in [9-11], which also have 
degenerate horizons that in the infrared source the 
flux. 

2. One may explicitly introduce bulk charged fields 
that source the electric field. If these charge carriers 



are fermionic then the eventual ground state of the 
system generally involves a gas of fermions filling 
(much of) the bulk geometry [12, 13] (see also [14- 
16]): this state was dubbed the "electron star" by 
[13]. If these charge carriers are bosonic then at 
sufficiently low temperatures the scalar generally 
condenses in the bulk, breaking the U{1) symme- 
try: this is the well-studied holographic superfluid 
[17], and will be discussed in the next section. 

It is of course possible to consider these two situations 
coexisting; i.e. there could be a horizon with a nonzero 
electric field, together with some charged matter outside 
the horizon that also contributes to the final boundary 
value of F^*. Indeed, at finite temperature this is generi- 
cally the case, and the recent construction of [18] realizes 
this explicitly at zero temperature. 




FIG. 1. Different ways to create a finite U(l) density in a 
field theory with a holographic description. Left: The electric 
field at infinity is sourced by a horizon. Right: The electric 
field is sourced by charged fields in the bulk. 

We will show that for a holographic system in the 
above setup, instead of (1.3) one has 

-L_Y.1^V. = P-^ (1.5) 

where the deficit A being given by the electric field flux 
through the horizon of the bulk geometry, i.e. given by 

A= (^-^V^F^\rh)^ . (1.6) 

Equation (1.5) means that the total Fermi surface volume 
only counts the charge density outside the horizon. 

This formula is not novel. When A = 0, it was pre- 
viously shown in specific examples in [22, 23]. In the 
form written above the formula was conjectured in [21] 
and also follows from the discussion in [19, 20, 24]. In 
particular an elegant and general proof for a holographic 
system with no horizon was given in [24]. Our deriva- 
tion is along similar lines as [24], but differs in some de- 
tails. We work in a more general situation, allowing for a 
horizon and the associated nontrivial widths of fermionic 
excitations. 

We now discuss the standard understanding of such 
a deficit, as articulated previously by various authors 
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[1, 19-21, 23, 25]. The gapless excitations which char- 
acterize a Fermi surface (as defined in (1.2)) are exci- 
tations of gauge invariant operators, and may be con- 
sidered "confined" degrees of freedom, corresponding to 
the quanta of mesonic bound states in the gauge the- 
ory. On the other hand, in holography [26], one thinks 
of the degrees of freedom associated with a horizon as 
being "deconfined" . Equation (1.5) then indicates that 
the deficit in the Luttinger theorem is associated with 
deconfined degrees of freedom. This appears to resonate 
well with the many-body results [4-6] mentioned earlier, 
where the Luttinger deficit is indeed due to an emergent 
gauge symmetry. Note that to make the identification be- 
tween a holographic system and these many-body states, 
one should identify the holographic gauge group with the 
emergent slave-particle gauge symmetry, in which case 
gauge- invariant fermionic operators can be interpreted as 
"electrons", and the deconfined horizon degrees of free- 
dom can be interpreted as "fractionalized" excitations. 
Defining 



A 



(2^) 



(1.7) 



it is thus tempting to speculate that in a general system 
(not restricted to those with a gravity dual) the deficit 
A provides a measure for the deconfined (or in the slave- 
particle context "fractionalized" ) degrees of freedom. For 
holographic systems A then has an elegant geometric de- 
scription (1.6) in terms of electric flux at the horizon. 



B. Deficit in the Magnus force 

We now consider the case where the U{1) symmetry 
is broken by a charged condensate, i.e. we study a su- 
perfluid phase. Before turning to holographic theories, 
consider creating a vortex excitation in a conventional 
superfluid and moving it slowly through a closed loop of 
area A. In this process the vortex will pick up a Berry 
phase 9. In a conventional superfluid this phase is pro- 
portional to the total charge enclosed by the loop, i.e. 



2iiA 



(1.8) 



This Berry phase actually means that there is a trans- 
verse force Ft on a moving superfluid vortex. This is a 
well-known (if somewhat controversial) force, and in the 
literature is called the Magnus (and/or lordanskii) force. 
In an ordinary superfluid one has 



27T 



compared with (1.8) and (1.9), which has a similar origin 
to the Luttinger count deficit mentioned in the previous 
subsection. 

Keeping this in mind, we now turn to the holographic 
superfluid [17], in which the U{1) symmetry is broken by 
the condensate of a charged scalar field in the bulk. In 
general there is also a horizon carrying a charge density. 
Interestingly, in such systems we also find that 



27rA 



ip-A) 



and 



27r • ■ 



(1.10) 
(1.11) 



where A is again given by (1.6). That is, the Magnus 
force also only picks up the charge outside the horizon. 
In parallel to the discussion of last subsection, defining 



A = p- 



qFr 
27rv ' 



(1.12) 



it is again tempting to speculate that A gives a mea- 
sure of the deconfined charged degrees of freedom in a 
superfluid phase. 

It is rather intriguing that in the holographic con- 
text, the deficit (1.7) for the Luttinger theorem and the 
deficit (1.12) for the Magnus force of a superfluid vortex 
have precisely the same bulk origin in term of the elec- 
tric flux (1.6) at the horizon. As we will see in the main 
text, there is a common thread in the derivation of both 
deficits: it is the bulk Gauss's law, which relates the bulk 
electric field to hulk charge densities: 



V -E{r) = pbuik{r) . 



(1.13) 



-e-^Vjp 



(1.9) 



It is the bulk electric field (evaluated at infinity) that 
determines the boundary theory charge density, but it is 
actually the bulk charges that contribute to the various 
observables that we compute. The horizon provides a 
mechanism for these two quantities to not agree, as will 
be more clear as we proceed. 

The plan of this paper is as follows: in Section II we 
specify the sort of action that we will be working with 
and fix notation. In Section III we prove the modified 
Luttinger theorem (1.5) and apply it to various systems 
in the holographic literature. In Section IV we turn to the 
bosonic superfluid case. We provide a brief review of the 
physics of vortices in holographic superfluids and then de- 
rive the form of the Berry phase and Magnus force, again 
discussing how the results apply to various systems in the 
literature. Finally Section V provides a brief discussion 
on the possible interpretations of these results. Various 
technical details are relegated to two appendices. 



where q is the charge of the condensate. 

However, it was found in [4] that in certain exotic su- 
perfluid phases with an emergent Z2 gauge symmetry, 
the Berry phase and the Magnus force exhibit a deficit 



II. SETUP 

We will be holographically studying strongly-coupled 
d-dimensional field theories with a conserved U{1) cur- 
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rent j^. We assume very little about this field theory, 
except that it has a limit in which it can be assumed to 
be dual to semiclassical gravity governed by the standard 
Einstein-Maxwell action 



So 



,[A,g]= J d^+h 



-^(2^ (^-2^) -4^^^ 

(2.1) 

Here F is the field strength tensor of a U{1) gauge field 
A that is dual to the conserved current j^. is a bulk 
gauge coupling, which is assumed to be small. In each 
of the following sections we will add either fermonic or 
bosonic charged matter to this action. Note that we do 
not actually assume any particular form for the gravita- 
tional part of the action except that it support solutions 
with horizons and the usual sort of conformal boundary 
at infinity. The bulk gauge coupling can also be taken 
to depend on radius and so be the expectation value of 
a bulk scalar field, in which case our results apply to 
dilatonic systems, such as those studied in [9, 18]. 

We assume that the bulk metric is diagonal, depend- 
ing only on the holographic direction r and preserving 
translational and rotational invariance. 

ds'^ = —gttdt^ + Qrrdr'^ + gijdx^dx^ (2.2) 

Throughout this work M, N will run over all the bulk 
directions and /i, v will run only over the field-theory 
spatial directions. While the fermionic portion of our 
discussion will apply to any dimension, in the scalar por- 
tion we will restrict do d = 3, so that the superfiuid 
vortices are pointlike excitations. In the more interesting 
cases we will assume that the space-time has some sort 
of a horizon at r = r/^; this may be a degenerate or a 
finite-temperature horizon. 

In general, our notation for spinors is that of [30]. 
Gamma matrices with an underlined index (e.g. V-) are 
in an ort honor mal frame. 



III. BULK FERMIONS AND LUTTINGER'S 
THEOREM 

In this section we discuss the analog of Luttinger's the- 
orem (1.5) as applied to holographic systems with charge 
carried by bulk fermions. Thus we add to the geometric 
part of the action (2.1) bulk fermions. 



,[il),A] = - j d'^^^x^itjjiT^VM - m)tlj (3.1) 



Here V contains the spin connection as well as a cou- 
pling to the background gauge field with charge q. This 
bulk fermion il){r^x) is dual to a charged fermionic op- 
erator which we will call 0{x). On many gravitational 
backgrounds [22-24, 27-30] the correlation functions of 
O exhibit Fermi surfaces, i.e. near a shell in momentum 
space they takes the form 



uo — vpik — kp) + 5](co', k) 



(3.2) 



The precise form of the self-energy S(co',/c) depends on 
the infrared geometry of the gravitational background in 
question; see [31] for a review. We will show that for 
each U{1) gauge symmetry in the bulk there exists a sin- 
gle constraint relating boundary theory charge densities 
with these Fermi surface singularities. A very similar re- 
sult was recently obtained in [24] for a system with no 
horizon; our derivation is slightly more general in that 
we allow for a horizon and treat carefully the analytic 
structure and nontrivial quasiparticle widths associated 
with such a horizon. 



A. Proof 

We would like to relate excitations of the bulk fermions 
to the bulk electric field. As we are dealing with fermions 
this will be an intrinsically quantum-mechanical contri- 
bution, and so we need to perform a one-loop calculation 
in the bulk. Note that the previous literature on the 
electron star [13] should also be formally considered as a 
"one-loop" calculation, but in a fiuid limit in which the 
scales determining the local electronic density of states 
were far higher than the curvature scale of the geome- 
try; thus it is justified to treat each individual fermion as 
being arbitrarily well-localized, and the gas of fermions 
is treated as a fiuid with the d + 1-dimensional equation 
of state that is appropriate to fiat space^. We will not 
make any such assumptions, and will treat the radial de- 
pendence of the fermion wavefunctions exactly. 

We will treat the gauge field and metric classically. We 
perform most of our calculation at zero temperature in 
Euclidean signature. The analytic continuation is 



t 



-IT 



: IUJE 



iS ■ 



(3.3) 



where Se is the Euclidean action. We extend our results 
to finite temperature at the end. 

To begin, we simplify the spinor action by rescaling 
the bulk field as 



(3.4) 



As is well-documented [30], this form eliminates the 
explicit appearance of the spin connection, and the 
Lorentzian spinor action in momentum space becomes 

S^ = - j dr^/g;;^ (2n)d-^ ^'^^^ ^' r)D^(cj, k; r) 

(3.5) 

where the bulk differential operator D is 

z) = i(r^a, + v^(-m + ir^if^v^)), (3.6) 



^ This equation of state is the standard one from statistical- 
mechanics textbooks; of course if one wants to obtain it from the 
field theory of free fermions it does involve a one-loop fermion 
determinant. 
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with Kt = (— — At{r)),Ki = ki. This operator D wih 
have eigenfunctions: 



(3.7) 



Here the are defined as usual by specifying bound- 
ary conditions on the spinor fields at both ends of the 
spacetime. The choice of boundary conditions on spinors 
at the AdS boundary is discussed in [32]. We will dis- 
cuss the boundary conditions at the horizon as we need 
them. Note that eigenfunctions with = satisfy the 
bulk wave equation as well as the boundary conditions, 
and so are normalizable modes and thus signify a singu- 
larity in the boundary-theory Green's function (i.e. the 
correlation function of the dual operator O). 



Mr] 



FIG. 2. Feynman diagram representing one- loop fermion den- 
sity. 



We now want to compute the effect of the fermion on 
the gauge field. Formally this is accomplished by inte- 
grating out the fermions, i.e. we work in Euclidean space 
and define a contribution to the effective action T[A] for 
the gauge field via 



exp(-r[A]) = J[d^exp{-SE[^P,A]) = {detD), (3.^ 



where the determinant involves all modes in Euclidean 
signature, i.e. D is evaluated with u purely on the imag- 
inary axis. Thus the boundary conditions at the horizon 
for the eigenfunctions (3.7) are now fixed: we should de- 
mand that the wavefunctions be regular at the Euclidean 
horizon. 

We now want to compute the contribution of this quan- 
tum effective action to the classical equation of motion 
for Ar. Adding T[A] to the geometric part of the action 
(2.1) and varying with respect to A^ we find the following 



equation of motion: 



■ ST[A] 
SArir) 



^ SD 



(3.9) 



This is simply Gauss's law (1.13) in the bulk. The right- 
hand side is the pointwise charge density carried by the 
fermions, where in the last equality we have used (3.8) to 
express it in terms of the operator D. Note that this is 
the Feynman diagram shown in Figure 2. We now need 
to evaluate this right-hand side. This can be done in 
several ways. In this main text we present a somewhat 
slick derivation using locality in the radial direction, and 
in Appendix A we present a more pedestrian treatment 
using the spinor bulk-to-bulk propagator. 

To begin, note that at fixed (Euclidean) frequency and 
momentum D is an operator acting on functions ^(r) 
that depend only on the radial direction. A basis for 
this function space is provided by the position eigenstates 
|r), where we can normalize to have the completeness 
relations 



(r|/) 



1 



-_S{r-r') (3.10) 



Note that we assume that the basis is complete using only 
r > which is true provided that our function space 
is Euclidean frequency wavefunctions that are smooth at 
the horizon. Now D is almost a local operator in r, and 
furthermore due to its gauge-invariance we can trade a 
derivative with respect to A^ for a derivative with respect 
io uje to find 



ri 



8D 



r2 



qiri 



dD 



dujf 



r2)S{r-r2) (3.11) 



Inserting this representation into (3.9) and performing 
some trivial manipulations we find the useful relation 



SArir) 



D- 



dD 



duJE 



(3.12) 



This equation is valid at all points in the bulk but is 
somewhat formal and unenlightening; however if we in- 
sert this again into (3.9) and integrate from the horizon 
to infinity then we find a trace over D-^duj^D, which we 
can rewrite as a sum over the eigenvalue spectrum of D 
(3.7). The final expression is then 



duE d"^ ^ d 
-^Q I 7^7T^ZwrT2^^1ogAn(cj£;) =/ 



(27r) (27r)^ 



dujE 



(3.13) 



r 



We have now removed the radial integral from the prob- 
lem; the manipulations that follow are essentially those 
from regular field theory, although we will repackage 



them in a way that allows for easy generalization to mul- 
tiple Fermi surfaces. 

As we now explain, (3.13) has something of a topo- 
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logical character. Note that the integration contour goes 
straight up the imaginary uj axis. Imagine that we could 
close the contour on the left. Morally speaking, the in- 
tegral is then counting the number of zeros of the eigen- 
value spectrum in the left-half plane. Each of these zeros 
corresponds to a Fermi surface, and so the integral over 
k will count the volume enclosed by the Fermi surfaces. 
This would have been a precise statement if the had 
been analytic functions of cj; however they actually con- 
tain singularities on the real uj axis, and so to make this 
intuition precise we will need to work slightly harder. 



AAA/VWW\A^<W\/VWWNA/W-> 



FIG. 3. Contour manipulations used in evaluation of (3.13). 
Squiggles indicate non-analyticities in \{uj) on real axis. 

Thus consider deforming the integration contour as 
shown in Figure 3 so that it wraps the negative real uj 
axis. Now recall that \{ijo) is defined by analytic con- 
tinuation from strictly imaginary uj] in the upper half 
plane the continuation to real uo continues to infalling 
boundary conditions, defining a function A^(cc;), whereas 
in the lower-half plane it continues to outgoing boundary 
conditions, defining a different complex function X^{uj) 
[32]. These conditions do not agree on the real uj axis, 
leading to non-analyticities in A(co') there: thus we must 
use X^{uo) for the top half of the integral and A*^(co') for 
the bottom half. These functions are not unrelated: one 



can show that for fermions with time reversal symme- 
try A^(a;) = -A^(cj)*. The integral (3.13) can now be 
written as 

I = 2,J^j_^-^a,i.), (3.14) 

where ak{uj) is defined to be 

ak{uj) = ^argA^(a;) . (3.15) 

n 

Now note that the integral essentially counts the winding 
in ak{uj) from the origin to — oo. As this is a topological 
expression, the full dependence of this expression on k is 
expected to arise from the endpoint at uj = 0.^ Note fur- 
ther that Xn{oJ = 0) changes sign as k is moved through 
each Fermi surface: thus as k is increased a/e(0) jumps by 
—IT each time a Fermi surface is crossed. Thus we may 
write the integrand as 



(27r)^ 



-N{k) 



(3.16) 



where if we denote the set of Fermi surfaces by {kp^}, 
then N{k) is a piecewise constant function: 



N{k) 



No 

No-1 
No -2 

' ' ' 1 





k < k^P 



kjp <C <C kpi 



k 



(3) 



(3.17) 



where Nq is a constant. In this expresion we assume 
that N{k oo) = 0; the necessity for this will be shown 
shortly. Denote by Vi the volume enclosed by the i-th 
Fermi surface. Using the expression (3.17), we find 



(27r)^ 



{NoVi + {No - 1){V2 - Vi) + {No - 2){Vs - F2) + • • •) 



(27r)« 



(3.18) 



r 



Somewhat magically, across each pair of terms the term 
involving A^o vanishes, and the expression above reduces 
to the sum of Fermi surface volumes. Putting this back 
into (A12) and using the AdS/CFT expression (1.4) for 
the boundary theory charge density, we find the desired 
result, i.e. 



(27r)^ 



A, 



(3.19) 



We assume the behavior at infinity is not affected by k. 



where p is the total boundary theory charge density and 
A is an anomalous term that measures the electric fiux 
through the horizon. 



A = 



9f 



(3.20) 



As claimed above, the sum of Fermi surface volumes mea- 
sures only the charge outside the horizon. 

Now note from (3.18) that if N{k ^ 00) 7^ 0, then the 
integral would receive a divergent contribution propor- 
tional to A^(oo) multiplying the infinite volume in /c-space 
that is outside the last Fermi surface. This seems like a 
clear UV pathology that should be avoided: by demand- 
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ing that N{oo) = we are essentially asserting that there 
is no "Fermi surface at infinity". While this seems quite 
reasonable, this involves UV properties of the spectrum 
and we do not know how to prove this to be the case in 
general. Note from (3.17) that this requires that A^o be 
equal to the total number of Fermi surfaces. 

Before moving on, we outline a generalization; consider 
heating the system up to a finite temperature. If the 
field theory is in a deconfined phase we will have a black 
hole horizon in the interior; all excitations will acquire 
T-dependent widths, and all of the Fermi surfaces will 
become somewhat blurry. Nevertheless, there is still an 
exact Gauss's law at each point in the bulk, and so some 
analog of the above expression should still exist in the 
black hole geometry. 

It is shown in the Appendix that instead of (3.18) one 
instead finds 

(3.21) 

where t/j is the digamma function and the sum over a is 
now a sum over all of the complex quasinormal modes 
of the system ^a{k). To understand the role of this 
digamma function, it is helpful to consider the case when 
the quasinormal modes are real, in which case one finds 

(3.22) 

which is simply the familiar expression from Fermi liquid 
theory. Thus we see that the digamma function has the 
effect of smearing out the sharp expression seen in (3.18), 
and should be viewed as a generalization of the normal 
Fermi-Dirac distribution to complex frequencies. This is 
in the spirit of [33, 34], which showed that such gener- 
alizations play an important role in one-loop corrections 
to the thermodynamics of black holes. 

B. Applications 

We now discuss how the formula (3.19) relates to var- 
ious systems in the literature; most of the models dis- 
cussed may be viewed as different solutions to the theory 
defined by the sum of (2.1) and (3.1). Note that the 
derivation does not appear to be sensitive to the width 
of the excitation (provided that there is a normalizable 
state at = 0) and so this modified Luttinger's theorem 
should apply to non-Fermi liquid excitations as well as 
to stable quasiparticles. 

1. Confining geometries 

We begin with the recent construction of [24] , in which 
the bulk spacetime is explicitly cutoff in the infrared to 
model a confining dual theory. In this case there is no 



horizon, and all the charge is clearly carried by fermions; 
^ = and the Luttinger count is saturated. Indeed the 
formula (3.19) was previously explicitly derived in that 
work. It is argued there that the only gapless excitations 
are those associated with the Fermi surfaces, and this 
system is thus clearly a Fermi liquid. 

2. Electron stars 

We now turn to the electron star geometry [13]. In 
that solution there is no infrared cutoff; the space-time 
exists for arbitrarily large proper distance in the infrared, 
and is filled with a gas of fermions that sources the elec- 
tric field. In a fiuid limit for the fermions one can self- 
consistently solve for the radial dependence of the geom- 
etry, gauge field, and thermodynamic variables charac- 
terizing the fermions. This was done in [12, 13], where it 
was found that in the infrared the geometry exhibits an 
emergent Lifshitz scaling, i.e. it takes the form 

ds^ = LU + ^ A = ed- (3.23 

This geometry is invariant under the scaling 

t ^ Xt x^X^x a ^ Xa, (3.24) 

where z is an emergent dynamical critical exponent that 
can be related to the mass and charge of the bulk 
fermions. Note that the form of the gauge field is fixed 
by this scale invariance (up to the overall constant Cd). 
There is some sort of a degenerate horizon as we take 
(J oo; as we take <j ^ eventually the geometry 
crosses over to an AdS^i+i in the far UV. 

However, we can compute the electric field fiux through 
this horizon; from (3.23) we find A to be 

Cd d-l 

(3.25) 

9f 

which thus vanishes as <j ^ oo. Thus for all finite z there 
is no fiux carried by the degenerate horizon. By (3.19) 
we see that the sum of all the boundary theory Fermi 
surfaces must make up the full charge density. This was 
verified explicitly in [22, 23], who found by direct com- 
putation a dense set of boundary theory Fermi surfaces 
on the background whose infrared geometry is (3.23) and 
checked that indeed the sum of their volumes is equal to 
the total boundary theory charge density. Those calcu- 
lations were performed in a WKB limit; now we see that 
the result must hold even away from this limit. 

Note that even though the geometry is very different 
in this case than the confining example above, from the 
point of view of the charged sector and the Luttinger 
count, there is little difference and it still appears that 
the charged degrees of freedom are "confined" 



^ Note that due to the warp factor a z before the spatial section 
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3. Reissner- Nordstrom black hole 

In the previous examples the traditional Luttinger 
count was satisfied, as A above was 0. We now turn to the 
opposite extreme, which is one of the workhorses of finite- 
density holography: the well-studied charged Reissner- 
Nordstrom-AdS black hole [8]. This gravity background 
is a solution to the pure Einstein-Maxwell part of the 
action (2.1), with no fermions sourcing the background. 
The infrared geometry can be understood as the z ^ oo 
limit of (3.23), i.e. AdS2 x 



-dt^ + da^ 



c^dx^ 



A 



e/- (3.26) 



This horizon carries a charge that appears independent 
of the fermion field content, i.e. we have 



A = 



9l 



r2 



(3.27) 



At the level of classical gravity, this is the full contribu- 
tion to the boundary theory charge density; if we were 
to obtain this bulk action from an actual stringy embed- 
ding, it would scale as 0{N'^). The fermions play no 
role in this analysis; thus in the extreme large N limit, 
this is the opposite limit to the one considered above; all 
of the charge is behind the horizon, and the traditional 
Luttinger count is thus maximally violated. 

Nevertheless, if one probes this charge density with 
a bulk fermionic field (3.1), one finds that the correla- 
tion functions of the boundary theory fermion operator 
contains a small (i.e. 0(1)) number of isolated Fermi 
surfaces [27-30]. These take the form 



uj — vpik — kp) + < 



(3.28) 



with V a nontrivial scaling exponent [30] that is related 
to the scaling dimension of the fermion operator in the 
AdS2 region. One expects that these modes should also 
contribute to the total charge density. This section may 
then be viewed as a calculation of this small correction, 
which from (3.19) can be seen to be of 0(1) rather than 
0(^)~0(iV2). 

Over most of the parameter space the system also of- 
ten contains a set of incoherent excitations dubbed the 
"oscillatory region" by [27, 30]. Their existence can be 
traced to the AdS2 scaling dimension going imaginary. 
This set of excitations is quite peculiar; for example at 
zero temperature and frequency one finds that for all k 
less than a threshold value k^: 



lmGR{uj = 0;k <ko) ^ 0, 



(3.29) 



of the metric (3.23), any bulk degrees of freedom with a nonzero 
boundary momentum will have a divergent local proper momen- 
tum approaching a ^ oo. Thus for gapless modes near a Fermi 
surface, the geometry (3.23) essentially behaves as a confining 
geometry with a "soft" wall. 



i.e. even at precisely zero frequency there is nonvanishing 
spectral weight. Related to this, for k < ko one finds an 
infinite number of poles in the retarded Green's function, 
geometrically spaced along a line in the lower-half plane: 



nUk)=^oexp(ie-'^) 



n e Z+, (3.30) 



where is a UV scale related to the chemical potential, 
Qq ^ jii and Xk is related to the imaginary part of the IR 
conformal dimension of the fermion operator. 

One should now include the contribution of these exci- 
tations to the one-loop calculation done above. It is not 
at all clear how to do this in a controlled manner, but to 
obtain a crude estimate we can imagine going to finite but 
very small temperature. Then the spectrum is cut off at 
frequencies ~ T, and one finds that the total number 
of poles is N ^ log (^). Each of these poles will con- 
tribute to the answer via the formula (3.21); it is difficult 
to estimate the contribution of each pole, but assuming 
it to be a finite number # we find for the contribution to 
the oscillatory region 



=(T)~#log(^) 



(3.31) 



Note the T ^ divergence. This should be added to 
the "anomalous" contribution from classical gravity A to 
obtain the full field theory charge density as in (3.19). 
It is clear that at exponentially low temperatures this 
quantum correction will outweigh the classical horizon 
contribution and this calculation will break down: what 
the system is telling us is that we need to take into ac- 
count the effect of the fermions on the classical electric 
field (and thus also on the geometry). 

This conclusion was first reached in [12] via an argu- 
ment that is essentially a WKB limit of the calculation 
performed above. Once we include this backreaction ef- 
fect, the final geometry is that of the electron star (3.23). 
In particular, this set of incoherent excitations is resolved 
into a dense (but discrete) set of Fermi surfaces - these 
excitations suck the charge out of the horizon and sat- 
urate the Luttinger count. As emphasized in [23], this 
is an example in which going to low energies one finds 
a smooth crossover from a "deconfined" phase in which 
charges are described by a horizon to a "confined" phase 
where charges live outside the horizon.^ 



4' Dilatonic systems 

We now turn to systems that have a nontrivial kinetic 
term for the bulk gauge field, i.e. where the gauge cou- 
pling is essentially taken to be a function of a dilaton 



Here by "confined" and "deconfined" we only refer to charge 
degrees of freedom. 
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that can run: 
S[A]- 



-i J d''+^x^gZ{4>)F^ . (3.32) 



Generally 0(r) will have a nontrivial profile in the in- 
frared, allowing for more interesting infrared geometries 
than those discussed above. Note first of all that our 
derivation above never assumed that the gauge coupling 
gp was constant in r, and so all of the discussion can 
be immediately taken over, except that the electric field 
flux now depends explicitly on ^, e.g. instead of (3.20) 
we have 



(3.33) 



and so on in all relevant formulas. 

We discuss first the dilatonic black holes studied in [9] . 
These are solutions to the action (3.32) together with the 
Einstein-Hilbert term and a kinetic term for ^, but no 
fermions. They have an infrared geometry whose metric 
is Lifshitz with finite z as in (3.23), but the dilaton (j) 
runs away in the infrared in such a way that the flux 
through the horizon A is finite in the infrared, rather 
than the vanishing result found in (3.25). Thus despite 
the very different infrared geometry from the Reissner- 
Nordstrom case, from the point of view of the anomalous 
Luttinger count the systems appear very similar; they 
also maximally violate the Luttinger count. 

Finally, we turn to the construction of "fractionaliza- 
tion of holographic Fermi surfaces" in [18], who study 
a gauge field action of the form (3.32) together with 
fermions. In this construction one finds that generically 
there is a nonzero charge both behind the horizon (i.e. in 
A^O) and in a fermion density outside. The fraction of 
charge behind the horizon can be tuned as a function of 
boundary theory couplings, allowing one to interpolate 
between a phase that saturates the traditional Luttinger 
count and one that maximally violates it. This is the first 
zero temperature example in which none of the terms in 
(3.19) is zero, and so is a nontrivial application of this 
formula. 



IV. VORTICES IN HOLOGRAPHIC 
SUPERFLUIDS 

In the previous section we studied fermions that carry 
charge in the bulk, and we explained how there exists 
a field-theoretical observable (the sum of the Fermi sur- 
face volumes) that measures only this bulk charge, and 
is insensitive to the electric field flux behind the horizon. 
In this section we seek a similar observable for the case 
when the bulk charge is carried by bosons. In this case 
generically the bulk boson condenses, breaking the U{1) 
symmetry, and we find ourselves in a superfluid state [17]. 
See [35-37] for reviews discussing the physics of such a 
phase. It is perhaps not immediately clear what the rel- 
evant observable should be in this case. 



A clue is given by [4], who study phases of mat- 
ter with an emergent Z2 gauge symmetry. In the con- 
densed matter context, such emergent gauge structures 
in fermionic systems are often associated with violations 
of Luttinger 's theorem [5, 6]. There also exist superfluid 
states with such an order: one of the results of [4] is that 
vortex excitations in such systems provide an interesting 
probe of the charge density. 

To understand this, consider first a generic (conven- 
tional) many-body system at zero temperature in two 
spatial dimensions, in which the ground state is a su- 
perfluid that breaks a [/(I) symmetry. Such a state will 
generically have gapped vortex excitations, around which 
the fluid circulation is quantized in multiples of ^ , where 
q is the charge of the boson that is condensed. 

Now consider moving this vortex in a closed loop en- 
closing an area A. It is a standard result [38, 39] that the 
wavefunction of the system will pick up a Berry phase 
that is proportional to the area enclosed, and is 



27r 

q 



pA 



(4.1) 



where p is the total number density, which at zero tem- 
perature is stored entirely in the condensate. In other 
words, in a conventional system the Berry phase accu- 
mulated by the vortex counts the total charge density 
enclosed by the loop. 

There is a macroscopic manifestation of this Berry 
phase. To understand this, it is helpful to consider the 
familiar case of a charged particle moving in a magnetic 
field 5 = V X a. If this particle is moved through a closed 
loop F, it accumulates a Berry phase that is 



1 (p i 



~ ^ ^ dx • a 



qBA 



(4.2) 



i.e. the phase counts the flux enclosed by the loop. In- 
deed, reversing the logic of (4.2), the presence of such a 
phase implies a term in the classical Lagrangian that is 
linear in the velocity, which means that a particle moving 
in a magnetic field with velocity v feels a Lorentz force 
Fl = qv X B. By analogy, we conclude that the result 
(4.1) suggests that a vortex moving with velocity v feels 
a transverse force proportional to the density. 



27r 

q 



(4.3) 



This is called the Magnus force. This transverse force on 
a region of circulating fluid exists in any hydrodynamic 
system and is not restricted to superfluids, and in every- 
day life is responsible for such important effects as the 
motion of a curve ball and the lift force on an airplane 
wing. 

Now consider heating the system up to a finite temper- 
ature. The Berry phase ceases to be well-defined; nev- 
ertheless one expects the transverse force on a vortex to 
remain a well-defined observable. The precise form of 
this force is now a matter of considerable controversy 
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[40]. Basically, at nonzero temperatures the density has 
both a normal component pn and a superfluid component 
ps . This situation has been well-studied in the context of 
Galilean superfiuids, where the normal component comes 
from thermally excited phonons. It is not clear whether 
or not (and how much of) this normal component con- 
tributes: see Chapter 2 of [41] for a recent review. We 
do not attempt to review the debate here, but we point 
out that for that if the fluid is at rest, then there exists 
a set of arguments that lead one to conclude 

27r 27r 
FriTf = —{ps + Pn) e^'vj = jpe^'v, (4.4) 

i.e. both the normal and superfluid components con- 
tribute^. There exist other arguments that lead to other 
results: we mention this result here because we believe 
that our holographic calculation can be interpreted as 
lending some support to it, although in a somewhat sub- 
tle way. 

For relativistic superfiuids without an underlying par- 
ticle description the situation appears more confusing 
still, as the separation of the total charge density into 
a normal and superfluid component is no longer obvious. 
While there does exist a prescription for performing such 
a separation (see e.g. [42-44] for a recent discussion or a 
brief review in Appendix B), it does not appear that the 
holographic calculation of the Magnus force is sensitive 
to the difference between the two, as will be more clear 
as we proceed. ^ 

In any case, these results do suggest that superfluid 
vortices form a probe of the charge density. In the re- 
mainder of this section we will compute both the Berry 
phase and the transverse force felt by a vortex in a 
holographic superfluid with a horizon. Just as in the 
fermionic example above, we will find that only the 
charge outside the horizon contributes to these observ- 
ables. At the end we will comment on how our results 
relate to the controversy. 

A. Background on holographic vortices 

In this section we will be working with a 2 + 1- 
dimensional field theory so that vortices are point like 



^ The part proportional to ps is called the Magnus force; the part 
proportional to pn is conventionally called the lordanskii force. 

^ Note also that in the Galilean superfluid the symmetry that is 
broken - number density - is actually entangled with spacetime 
symmetry generators in a nontrivial way: e.g. the momentum is 
proportional to the number current, which is not the case in the 
relativistic system under study here. We are thus not entirely 
certain how easily results from the Galilean superfluid should 
be taken over to the relativistic setting, and one should perhaps 
use caution when comparing the two. Unfortunately, we are 
not aware of a careful relativistic fleld-theoretical Magnus force 
calculation to which we can compare our holographic results. We 
thank K. Balasubramanian for discussions in this regard. 



excitations. We will work with a simple bulk action, in 
which we add to (2.1) a bulk scalar field with charge 

S4^, M = - j d^x^g {{D^yDc^ + V{^U)) (4.5) 

This is the action of the well-studied holographic super- 
fluid [17], and it is well-known that at sufficiently low 
temperatures (compared to the chemical potential) the 
system generally orders into a superfluid state where the 
scalar (j) obtains a nontrivial bulk profile: 

^{r)=vo{r), (4.6) 

breaking the U{1) symmetry in the bulk. 

At{r) also has a nontrivial profile, indicating (via (1.4)) 
that the field theory state is at nonzero charge density. 
Some fraction of this charge is stored in the condensate 
outside the horizon. At finite temperatures, there is also 
a nonzero electric field flux through the horizon. These 
both contribute to the total charge density. 

The detailed thermodynamics at very low tempera- 
tures depends on the bulk potential, but we will remain 
at finite temperature and assume that the potential is 
such that this condensation occurs. All of our calcula- 
tions will be done in a probe limit where the scalar and 
gauge field do not backreact on the geometry; however we 
do expect the essential features of our discussion to sur- 
vive the inclusion of back reaction. One may assume the 
background geometry to be just the AdS-Schwarzschild 
black brane; we have a horizon at radius r = r^. 

Now on general grounds whenever we break such a 
U{1) gauge symmetry we expect to find a state corre- 
sponding to an Abrikosov string where the phase of the 
condensate winds through 27r around a line in the bulk. 
In holographic models such solutions have been explicitly 
constructed [45-48] by numerically solving the relevant 
partial differential equations. Many features of these so- 
lutions may be understood analytically: imagine such a 
vortex line localized in the field theory directions at a 
location xq, stretching from the AdS boundary down to 
the horizon at r = r/,,. Let us use polar coordinates (p, 0) 
for the field-theory spatial directions, centering the origin 
at X = xq. We expect the bulk fields to take the form 

0(x) = v(r,p)e^"(^) Ae(p^oo)^-^o(-] (4.7) 

Q \pj 

The profile v{r^ p) will be fixed by bulk dynamics; it must 
vanish at the vortex core p = and far from the vortex it 
should approach the background profile (4.6). As usual, 
the condition on Aq is found by appealing to finiteness 
of the energy far from the vortex by demanding that the 
integral of the term g^^{dea — qAq)'^ in the energy of the 
vortex be finite (see e.g. [49]). This means that the flux 
enclosed in the string is 

$ = (fds-A = ^ . (4.8) 
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Note that outside the vortex core nothing is sourcing Aq , 
so it is essentiahy "pure gauge", where the quotes indi- 
cate that the gauge parameter of which it is a derivative 
is not a single valued function. We can rewrite it as 

Am{x) = ^OmMx^xq) A{x;xo) = arg(f-fo) (4.9) 

Thus as we move in a circle around the vortex A sweeps 
through 27r. 

It is clear that the field-theory interpretation of this 
configuration is a pointlike vortex excitation located at 
xq. There are some subtleties here associated with 
boundary conditions, which we review in Appendix B. 



B. Berry phase 

We now attempt a naive gravity computation of the 
Berry phase accumulated by this vortex if we move it 
very slowly in a closed loop, i.e. make xo{t) a slow func- 
tion of t, as in Figure 4. The final answer is simple to 
understand; the vortex carries magnetic fiux If we 
move it through a closed loop in the bulk, all the bulk 
charge density that is enclosed in the loop will encircle the 
fiux tube and contribute a phase via the usual Aharonov- 
Bohm effect. Thus we naively expect the Berry phase to 
only notice the charge outside the horizon. 



limit of adiabatic transport the value of JT* is 0, and so 
from (4.9) we have 



dtj^A{x;xo{t))J* (4.11) 



Now consider this expression carefully. The integral over 
t simply measures the winding of A as we move the vortex 
in a closed loop. If x is inside the closed loop xo(t), then 
A winds through 27r. If it is outside the closed loop then 
the winding is 0. Thus the integral receives contributions 
only from points inside the loop, and is 



27rA 



/ 



(4.12) 



where A is the transverse area enclosed by the loop, i.e. 
the phase is proportional to the bulk charge enclosed. 

We now relate this to a boundary theory quantity. 
Consider the bulk Gauss's law: 



dr 



9l 



-gJ' 



(4.13) 



Again, integrate this equation from the horizon to the 
boundary and use (1.4); we find then that the Berry 
phase is 



— Si^ 



27iA 



(4.14) 



FIG. 4. Moving vortex slowly in a closed loop. The Berry 
phase counts the bulk charge enclosed by the loop in the bulk. 

We will compute this phase by evaluating the bulk 
Lorentzian action on the adiabatic trajectory. First, we 
note that the key term in the action is 



int 



-I 



d'^x^/^AuJ 



M 



(4.10) 



Here is the hulk current, built from the usual Noether 
procedure applied to (4.5). Now there is a background 

value of A^^^ that is sourced by the superfiuid profile; we 
will neglect this term as it does not change as we move 
the vortex, and from now on one should assume that At is 
purely the vortex contribution (4.9). On the other hand 
the background value of will be important. In the 



where as before A is the electric fiux through the hori- 
zon (3.20). This is the claimed result. Note that this is 
precisely the same pattern that we find in the Luttinger 
count violation in the fermionic case: the bulk observable 
is sensitive only to the charge outside the horizon, and so 
A should be viewed as an anomalous contribution. 

Let us now turn a critical eye to this computation. 
From the bulk point of view, we assumed that the string 
was moving rigidly as we moved it from the boundary. 
However in a spacetime with a horizon, this is clearly im- 
possible. No matter how slowly we may move the string 
at infinity, from the point of view of an observer at in- 
finity, the movement will never propagate through to the 
horizon; the pulses of movement will slow down and ap- 
pear to freeze just outside. There is some confusion as to 
how we should deal with the interaction of the string with 
the horizon. This confusion is dual to the field-theoretical 
statement that in a system with a horizon and the asso- 
ciated gapless degrees of freedom, there is no such thing 
as "adiabatic" transport, and so the idea of the Berry 
phase itself is not very well-defined. 

Thus while the calculation illustrated above is enter- 
taining, it should not be taken too seriously. 



C. Magnus force 

We will now compute an observable that should be 
well-defined even at finite temperature; the transverse 
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FIG. 5. Sending ripples down vortex string to compute Mag- 
nus force in linear response. 



force on a vortex moving with a velocity v. We will do 
this by sending low frequency ripples down the vortex 
string into the horizon as in Figure 6 and computing the 



force that must be supplied to sustain such a ripple. We 
shall see that this can be formulated as a linear response 
problem. To do this we need to describe dynamics along 
the vortex, i.e. we seek an action that describes ripples 
in the collective coordinates characterizing the solution 
(4.7). We assume that such an action is given by 



= Sng + Si 



int 



(4.15) 



where Sng is the Nambu-Goto action 



S=-J Tp{X)^det {daX^dBXN)gMN{X) 

(4.16) 

Here we treat the core of the vortex as a stringlike ex- 
citation with a tension Tp that can vary along the ra- 
dial (holographic) direction. We pick coordinates on the 
world sheet (r, t) that correspond to the embedding co- 
ordinates in spacetime; then the quadratic term for fluc- 
tuations in the transverse X* is 



quad 



(4.17) 



Let us now consider the boundary conditions on the 
string. Consider the canonical momentum of the string 
coordinate with respect to a foliation in the r-direction, 
i.e. at the quadratic level we have 



(4.18) 



Note that if we are applying a force Fi on the string at 
the UV boundary, the value of this canonical momentum 
at infinity should be equal to the force applied, i.e. 



n,(r 



(4.19) 



I 

We turn now to the other term in the action. As be- 
fore, Sint describes the interaction of the vortex with the 
background charge density and is given by (4.12), except 
that now Am depends explicitly on the location of the 
string, i.e. 



A 



1 d 



M 



- A(x; X) A(x; X) = arg - 



(4.20) 

Note now that as only J'* is nonzero, the action is only 
sensitive to A^, which we write via the chain rule as 



When Fi = this is the normal Neumann boundary con- 
dition. This forced result is familiar in the context of the 
classic string dragging calculations in AdS/CFT [50, 51]. 
Thus to compute the force on the vortex we should com- 
pute this canonical momentum at infinity. 

I 



Thus the term in the action may be written as 



(4.21) 



Sir) 



I 



drdty/grrgttX^cii 



-I 

q J 



d 

■A{x',X)j' 



(4.22) 



Here we have defined an effective vector potential a^, which couples to each bit of the string the same way as a vector 
potential couples to a charged particle. Note that this action appears nonlocal, as it appears to involve an integral 
over the full field-theory spatial directions, and not just an integral over the string worldsheet''. We will show that 
the equations of motion are however local. 



Note also that we are also neglecting any effects from the core 
of the vortex; one should assume that the fluctuations X* are 



larger than the core but still small enough to be characterized 
by linear response. 
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Now we vary Sint- We find after some algebra 

This equation is valid for any vector potential (note that it is simply taking the curl). In our case, appears to 
be the derivative of a scalar function A(x; X) and naively the curl would vanish; but of course the scalar function is 
not single- valued, and now we use the following identity: 

[Sx, Sy]A(x; X) = (27r)(5(2) - (4.24) 

to convert the integrand into a delta function in the spatial directions. Thus the equations of motion are local on the 
string worldsheet. The final equations of motion are 

Ott 

dr {Tp^g^-drX') + Tp^g^'d^X' + e^'^V^j'dtX^ = . (4.25) 



^SXY - 5YX) [^-^ 



(4.23) 



This is the massless scalar wave equation plus an extra 
mixing term. 

We may decouple these equations by introducing the 
linear combinations and their canonical momenta 



X±=X±iY 



after which we find 



n± 



^S,X±), (4.26) 



drU± 



27r 



gg''cj%X^ ± cj^^j'X^ = (4.27) 



The above manipulations were valid on any background. 
We now specialize to a finite-temperature horizon. An ef- 
ficient method for the solution of equations such as (4.27) 
on such backgrounds was described in [52]. We do not 
review the method here; we just point that out a conve- 
nient object to consider is the following ratio: 



iujX-[ 



(4.28) 



It is shown in [52] that horizon boundary conditions force 
this to be a constant (call it cr{rh)) at the horizon. In our 
context, when evaluated at infinity it will relate the ap- 
plied force to the displacement of the vortex via (4.19)^. 
From (4.27) we find the radial evolution equation: 



drcr± — ± 



27ri 




(4.29) 



and now finally we go to the low- frequency limit, which 
lets us discard the second term. We can then immediately 
integrate to find the following answer 



Cr±(oo) = g{tu) ± 



27ri 



a{rh)±iM . 



(4.30) 

M counts the charge outside the horizon; as usual, we 
use the bulk Gauss's law (4.13) to express it in terms of 
field theory quantities as 



M 




27r 



(p-A). (4.31) 



The interpretation of this equation is far more clear if we 
convert back to the (X, F) basis and use (4.19) to relate 
the canonical momentum at infinity to the applied force. 
We then find 



.•.•27r 



(4.32) 



We see that this off-diagonal component of the force is 
proportional to the velocity, and is precisely the Magnus 
force. There is also a diagonal component proportional 
to cr(r^) that is purely dissipative and that we have not 
written down. (4.32) is the desired result. 



D. Applications 

We now discuss how these results apply to various sys- 
tems in the literature. 



If X* were a bulk gauge or gravitational mode, this same object 
a{r — oo) would have been a transport coefficient [52]. Indeed 
there is a great deal of formal similarity between the Magnus 
force computation described here and a holographic computation 
of a Hall conductivity; it would be interesting to understand 
if this is related to some form of particle-vortex duality in the 
boundary theory. 



1. Finite temperature superfluids 

We first study the usual holographic superfiuid, i.e. the 
model studied in [17]. For reasons to be clear below we 
call this the deconfined holographic superfiuid. All of the 
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results above apply directly to this model; we see that the 
coefficient of the Magnus force (4.32) in fact agrees with 
the naive Berry phase calculation (4.14). The coefficient 
of this force thus directly measures the charge outside the 
horizon. At finite temperatures this will disagree with the 
total charge density by the anomalous term A; we point 
out that the structure of the violation is identical to that 
seen in the rather different Luttinger count calculation 
for fermionic systems. 

We now relate this to the controversy discussed earlier. 
As mentioned earlier, in conventional Galilean superflu- 
ids it is a somewhat controversial matter whether or not 
the normal component of the number density contributes 
to the transverse force [40]. The normal density in ques- 
tion is a thermally excited phonon gas, which at finite 
temperatures coexists with the superfluid density that is 
locked in the condensate. It is interesting to see how our 
holographic computation relates to this matter. 

First, note that in holography it does not seem that 
the distinction lies in whether or not the charge density 
is "normal" or "superfluid" , but rather in whether it is in- 
side or outside the horizon, i.e. whether it corresponds to 
gauge- invariant or deconfined degrees of freedom. These 
need not be correlated; one could imagine putting some 
sort of extra charged bulk quanta on top of the super- 
fluid background, thus giving us a "normal" component 
that is nevertheless outside the horizon. Such an extra 
contribution would be a gauge- invariant density carried 
by a charged collective mode of the system, and so would 
be somewhat analogous to the phonon gas that is studied 
in conventional Galilean superfluids. This would clearly 
contribute to the expression (4.32); thus in holographic 
systems we believe one can say cleanly that any density 
- "normal" or "superfluid" - contributes in exactly the 
same way, provided it is outside the horizon. Indeed there 
is no need to separate them in the calculation. 

Of course in our holographic systems there is a new 
component to the charge density that does not con- 
tribute: the charge carried by the horizon. This should 
be roughly thought of as being carried by N'^ deconfined 
gauge-charged degrees of freedom in the dual field the- 
ory. There are no analogs of such degrees of freedom in 
the conventional superfluid; thus in any comparison with 
conventional systems the term A should be set to 0, and 
our holographic computation can then be viewed as sup- 
port for the many-body calculations leading to (4.4), in 
which both the normal and superfluid components con- 
tribute. 

We stress that one lesson here is that the deconfined 
holographic superfluid is thus not a conventional super- 
fluid. This point has been made before: at zero temper- 
ature the superfluid geometry is generally Lifshitz, with 
extra gapless modes that have no place in a traditional 
superfluid. At finite temperature one is essentially excit- 
ing these gapless modes, and we now see that one sharp 
observable to which they contribute is the coefficient of 
the Magnus force. 



2. Confining superfluids 

Conversely, if one can gap out the gauge theory dynam- 
ics (for example by considering it in a confined phase), 
then one expects to find a conventional superfluid. This 
problem has been studied in [53, 54], where the confined 
geometry is given by the AdS soliton [26]. 

We briefly discuss the construction below. Consider 
taking a (3+l)-dimensional gauge theory and compacti- 
fying one spatial direction on an parametrized by an 
angle ip. At low energies we are now working with a 
(2+1) dimensional gauge theory, and depending on the 
boundary conditions of fields around the we may ex- 
pect the theory to confine. A/" = 4 Super Yang-Mills can 
provide a precise realization of this [26], but we will be 
working in a bulk gravitational description with extra 
fields added, and so we will not have a precise CFT dual. 
The confinement is realized geometrically in that the 
direction shrinks smoothly to zero size in the interior, 
ending the geometry at some radial coordinate r = r*. 
This geometry is the AdS soliton. 

We may now consider adding to this bulk gauge fields 
and charged scalars to realize a confining superfluid; the 
phase diagram was mapped out in [53, 54], and there ex- 
ists a confining superfluid phase extending down to zero 
temperature over a range of parameters. This has no ex- 
tra gapless degrees of freedom and so should be dual to a 
conventional superfluid. It is instructive to see how our 
arguments above apply to this phase. 

First, it is necessary to identify the "vortex". The UV 
theory is (3+1) dimensional, and so vortex excitations are 
strings and not points. To obtain a point like excitation 
in the effective low-energy (2 + 1) dimensional theory, we 
should arrange for the line to wrap the compact S^. In 
the bulk the shrinks to zero size; thus the vortex sheet 
in the (4+1) dimensional bulk should be thought of as a 
membrane filling the entire (r, tjj) plane at a fixed value of 
the two remaining spatial coordinates xq. Now we have 
no horizon, and the Berry phase arguments above may 
be made precise, with ^ = 0. 

Similarly, we can repeat the linear response calcu- 
lation, except with different boundary conditions. At 
r = r* we are at the origin of polar coordinates, and 
so all radial derivatives of bulk fields should be 0. Thus 
from (4.28) we see that the new boundary condition is 
simply cr±(r*) = 0. Repeating the calculation we find 
the same result (4.32), except that now there is no dissi- 
pative force at all, and ^ is 0. This is the conventional 
superfluid result. 



3. Speculations at zero temperature 

We now return to the deconfined case. Here it would 
be very interesting to extend our calculations to zero tem- 
perature. We are hindered by the fact that as of yet there 
has been no explicit construction of the vortex soliton at 
zero temperature. This may be nontrivial, as generically 
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at zero temperatures one finds Lifshitz solutions; tiiese 
contain mild "singularities" [55-57] which lead local bulk 
observers to feel diverging tidal forces in the far infrared. 
It is not obvious to us that the vortex string will main- 
tain its coherence arbitrarily far in the infrared, and the 
fate of the vortex state thus seems to be an interesting 
topic for future study. 

If we assume that the soliton remains coherent, then 
we find the interesting result that at zero temperatures 
A drops to zero by the arguments surrounding (3.25). 
This is then somewhat similar to the situation seen in 
the electron star, in which we see no anomalous Lut- 
tinger count. As pointed out in [18], however, if we 
study systems with dilatonic couplings such as (3.32) 
except with a charged scalar condensate (rather than 
charged fermions), it seems that there might exist su- 
perfiuid analogs of the states studied there. In this case 
some charge remains behind the horizon and coexists 
with some charge outside stored in the scalar condensate. 
This system would then exhibit an anomalous Magnus 
force even at zero temperature. 

V. CONCLUSION 

In this paper we studied field-theoretical observables 
that should be thought of as probes of the charge density 
of the system. For fermionic systems, the observable is 
the sum of the volumes enclosed by all visible Fermi sur- 
faces. For bosonic systems, the probe is the coefficient 
of the Magnus force felt by a superffuid vortex. In both 
cases in "conventional" many-body systems, this probe 
is sensitive to the full density of the system, as is formal- 
ized for fermionic systems by the traditional Luttinger 
theorem (1.3). In holographic systems, however, we find 
a different result: for fermionic systems we have a modi- 
fied Luttinger theorem 

^^Y.V, = p-A, (5.1) 

and for bosonic systems we find 
. • 27r 

F'=e'^^(p-Ay (5.2) 
Q 

where p is the net U{1) charge density, and A is an 
anomalous term that in the holographic description is 
the electric field fiux sourced by the horizon. This term 
is not present in conventional systems. Fascinatingly, the 
horizon provides a mechanism for the violation of the 
traditional Luttinger theorem and its superfiuid analog. 
Furthermore, the violation takes exactly the same form 
in both the bosonic and fermionic case. 

We now note several facts that we find interesting re- 
garding this result. First, from the gravity perspective, 
we now have a field-theoretical definition of the charge 
carried by the horizon. Given a field theory, simply mea- 
sure the right-hand side of (5.1) or (5.2) and subtract it 



from the total charge density p to determine A. Note 
that this works at finite temperature, even in the spinor 
case when there are no well-defined Fermi surfaces; the 
left-hand side is just given by the somewhat messier ex- 
pression (3.21), which can in principle be found given a 
perfect knowledge of the boundary theory spinor Green's 
function. 

Before further interpreting these results, we should dis- 
cuss the possibility that the term A may be an artifact of 
the gravitational description (or equivalently, an artifact 
of the large N approximation we are using) , and that the 
true field theory does not have such an anomalous term. 
It might happen that at finite coupling or N the charged 
horizons discussed above would be replaced by some ob- 
ject in the gravitational description that would saturate 
the Luttinger count or contribute to the Magnus force. 
At this point we should discuss finite-temperature and 
zero-temperature horizons separately. If we are at finite 
temperature then the spacetime is completely smooth 
and we do not expect any such effects to change the 
results presented here. On the other hand, a zero tem- 
perature horizon is usually associated with an infrared 
subtlety such as a ground-state degeneracy or an IR sin- 
gularity, and here it seems quite possible that finite- 
or finite-A effects could significantly modify the physics. 
Earlier in Sec. IIIB3 we saw such an example: when 
the AdS2 dimension of a fermion becomes complex, the 
horizon of an extremal Reissner-Nordstrom black hole is 
replaced at finite A^ by an "electron star" which saturates 
the Luttinger count. 

While this is certainly a logical possibility, from now 
on we will assume that the effect is real and that an 
anomalous term A truly exists in the dual field theory. 
In that case we should attempt to make contact with the 
many-body systems that are known to exhibit a similar 
violation [4-6], in which one finds anomalous terms of 
precisely the form (5.1) and (5.2) in both the Luttinger 
count and the Magnus force. ^ As mentioned earlier, the 
existence of the anomalous terms is related to the exis- 
tence of a deconfined gauge symmetry. 

In fact, in the discussion of [4] of fractionalized Fermi 
liquid and superfiuid phases, there is a common thread 
in the derivation of the deficits for the Luttinger count 
and for the vortex Berry phase. In those derivations one 
finds relations similar to (5.1) and (5.2) by considering 
the system on a torus and threading 27r fiux around one of 
the cycles while tracking the flow of momentum. Essen- 
tially the conventional charged degrees of freedom (e.g. 
fermionis quasiparticles or superfluid vortices) pick up 
some momentum under this threading of flux, and if this 
was the whole story one would flnd the conventional Lut- 
tinger count or Berry phase. However in a fractionalized 
phase topological excitations of an emergent deconflned 



^ Other examples include those in [1], where the deficit is carried 
by some gauge variant spinon Fermi surfaces, which cannot be 
probed using gauge invariant operators. 
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gauge theory can also carry some of this momentum, re- 
sulting in an anomalous term similar to A. 

Fascinatingly, in the holographic considerations of this 
paper, there was also a common thread in the derivations 
of both results: the observables that we compute are sen- 
sitive to bulk charges^ while the actual field theoretical 
charge density is sensitive instead to the boundary value 
of the bulk electric field. The relation between these two 
quantities is provided by Gauss's law: 



V • E{r) = phuikir) . 



(5.3) 



However, if we have electric fiux emanating from behind 
a horizon, then there is a component of the boundary the- 
ory charge density that has nothing to do with any bulk 
charges, resulting in a nontrivial A that is also clearly re- 
lated to the existence of deconfined gauge fields. We thus 
speculate that the same mechanism is at play both in the 
holographic and condensed matter many-body systems, 
even if the respective derivations take rather different 
routes, and in both cases A provides a measure of the 
deconfined charged degrees of freedom. It is important 
to fiesh out this connection further. 

We note that it may be possible that the horizon charge 
degrees of freedom can be probed by some other observ- 
ables. Recently it has been proposed in [58, 59] that the 
horizon charge degrees of freedom in certain Einstein- 
Maxwell-dilaton systems may refiect "hidden" Fermi sur- 
faces which can be probed via entanglement entropy, but 
not through the correlation functions of gauge invariant 
operators. 

Further study will be required to understand precisely 
how such holographic finite-density phases fit into our 
conventional understanding of quantum matter. One 
might hope that such study will shed light both on the 
physics of horizons and of possible phases of matter. 
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Appendix A: Technical details of spinors 

1. Green's function evaluation of charge density 

Here we perform an alternative derivation of the one- 
loop charge density that uses an explicit Green's function 
representation for the spinor density of states. The final 
result is the same as (3.13) in the main text. We first 



review some technical aspects of spinor bulk-bulk propa- 
gators. Similar considerations can be found in e.g. [60]. 
The Lorentzian spinor action is 



giiPir^VM-m)^P 



(Al) 



By rescaling the spinor the action can be written in mo- 
mentum space as 



dk 



(A2) 



27r (27r) 

where the bulk differential operator D is 

D = i {T^dr + (-m + iT^K^^^)) 



(A3) 



with Kf = {—u — At{r))^Ki = ki. To define a propaga- 
tor we will need to invert this operator D. D will have 
eigenfunctions: 



D^n{r,uj,k) = An^n(r,cj,A:) 



(A4) 



To make this equation well-defined we need to specify 
boundary conditions on the fields both at the horizon 
and at infinity. At infinity we demand that the spinor 
fields be normalizable in the usual AdS/CFT sense, as 
discussed in e.g. [32]. To understand the horizon, let us 
begin in Euclidean signature, i.e. uj is entirely imaginary. 
In this case demanding that the equation be regular at 
the horizon forces the solution to die away there and 
completely fixes the boundary conditions. 

We now define an inner product between functions as 



/ 



dr <^T^^ 



(A5) 



D is not quite hermitian under this inner product; in- 
stead, one may show that at complex frequencies uj we 
have instead 



(A6) 



Nevertheless, this condition is enough to prove the or- 
thogonality of the eigenfunctions (A4), 



/ 



dr<S^{LU*;r)T^^„{oj;r) = S„ 



(A7) 



We are not certain how to prove completeness of the basis 
(A4), but if we assume it then the relevant completeness 
relation is 

J2 rWnito*; r')r^- = d{r - r') . (A8) 



One can check consistency by integrating both sides of 
this relation against another eigenfunction and using 
(A7). Given these relations, it is easy to see that the 
bulk-bulk Green's function is 
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By acting on this expression with D and using (A8) it is 
simple to verify that G satisfies 

DG{uj, k] r, r') = 6{r - (AlO) 

which is its defining property. We stress again that G{uj) 
has been defined by analytic continuation away from Eu- 



clidean cj, and so in the upper-half plane it will coincide 
with the retarded propagator, while in the lower half- 
plane it will coincide with the advanced propagator. 

We now use these expressions to rederive (3.13) using 
a Green's function representation. 



Our starting point is the variation of (3.8), which yields 



SAr{r) 



= -TrD" 



5D 



= +iq 



I 



dUE d^~^k 



27T (27r 



\d-l 



r^^/97r tr(r-G(ia;£;, k\ r, r)). 



(All) 



where the lowercase trace is over spinor indices. Note that this is essentially the one-loop Feynman diagram shown in 
Figure 2. This equation, while valid at arbitrary points in the bulk, is somewhat unenlightening; we now insert this 
into (3.9) and then integrate both sides from the horizon to infinity to find 



9f 



dr 



dujE d^ ^k 
^(27r)^-i 



^9' 



fg^ ir(V-G{iujE, k; r, r)) = / 



(A12) 



We now need to perform the radial integral to determine / explicitly. We first insert the eigenvalue decomposition 
(A9) to write the integrand as 



dr 



dujE d^-^k 
(27r)^-i 



■^n{iiOE;r)-^u{-iiOE;r)T^- 
A„(w) 



(A13) 



We now note that the coupling between the spinor and the gauge field takes the special form ^!{iuo — qAt)^ . Thus 
one can essentially trade couplings to the gauge fields for derivatives of the bulk wavefunctions with respect to uo] this 
manipulation is closely related to the Feynman-Hellman theorem. More precisely, we take a derivative with respect 
to uje of the defining equation (3.7): 



(A14) 



Now the first term in this expression is 



(A15) 



which is precisely the coupling leading to the F- in (A13). Thus we may insert (A14) into (A13) and then integrate 
by parts to transfer the D — \n io the other wavefunction factor, annihilating it. We then find only 



I = iq 



dujE d^ ^k ^dujE^n 
^(27r)^-i^ An 



dr'^ni-ioJE] r)T-'^n{ioJE, r) 



(A16) 



The radial integral is now done by appealing to the or- 
thogonality relation (A7). We find finally 

which is (3.13) in the text. 



2. Finite temperature Luttinger theorem 

We now present the manipulations involved in deriving 
a finite-temperature "Luttinger theorem", i.e. an expres- 



sion relating the charge density outside a finite tempera- 
ture horizon to a property of the spinor Green's function. 
We begin with the expression (3.13), except at finite tem- 
perature, so that the integral over frequencies becomes a 
Matsubara sum: 

where the fermionic Matsubara frequencies are 

uon = (2n + l)7rT n G Z (A19) 

We now perform some standard finite-temperature field- 
theoretical manipulations. We first write the Matsubara 
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sum as a contour integral: 



(A20) 



Now 



we use the identity 7rtanh(7rx) = 
ix) — i/j (^^ — ix)) to rewrite the hyper- 



bohc tangent and get the fohowing representation of the 
Fermi distribution: 



where the function f{oj) has poles at ah the Matsubara 
frequencies, and so is 
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ex.p{/3uj) + 1 
The relevant contour is shown in Figure 



f (tanh (^^) - 1 ) (A21) 



FIG. 6. Integration contour used in (A20). 



Putting this in we find the formula 



2tt J {2ny 



(A22) 



1 



2 \7r V V2 27T J ^ \2 27r 

(A23) 

This is useful because the digamma function ip only has 
poles at negative integer values of its argument; thus the 
first digamma has poles only in the upper-half plane, and 
the second has them only in the lower- half-plane. So we 
now do each digamma integral separately. 

For the first one, we deform the contour downwards 
into the lower-half plane, as in Figure 7. We are now 
working with the analytic continuation of Xn{^) to the 
lower-half plane; this has a collection of zeros at complex 
frequencies uj — Qa{k)^ each of which is a quasinormal 
mode. The integral over u will pick up a contribution 
times (— 27rz) from each of these zeros^^. We do the same 
thing for the other digamma, except that we now de- 
form the contour upwards, picking up a contribution from 

We find once the dust settles 



2 2n 

(A24) 

This is (3.21) in the main text: its interpretation is dis- 
cussed there. 



\UJ 




Q 



FIG. 7. Contour manipulations discussed above, for digamma function with poles in the upper-half-plane. Gray circles indicate 
zeros of X^{uj), which are present if X{(jj) is analytically continued into the lower- half-plane. Note a similar manipulation is 
done on the other digamma function. 



The minus sign is because the orientation of the integral is back- 



wards in the lower-half plane 
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Appendix B: Technical details of vortices 

Here we review some details about boundary condi- 
tions on the vortex soliton that were glossed over in the 
main text; some of these considerations can also be found 
in [47]. Recall that the bulk action is 



S 



I 



-9 I -^F' 



{D<j>)^D<j> + V {<!>'' <l>) 



(Bl) 

We will use r for the bulk holographic direction and will 
parametrize the field theory spatial directions by polar 
coordinates (p, 0). 

In the superfiuid phase the condensate (j){r) has a pro- 
file; it is large in the interior and dies away at infinity so 
that it represents a normalizable state in the theory. In 
the bulk the symmetry is gauged by A; thus there is some 
sense where we have a superconductor in the deep inte- 
rior that is crossing over into an unHiggsed phase near 
the boundary. If we now want to construct a vortex on 
top of this background superfiuid, we put an Abrikosov- 
Nielsen-Olesen string on top of this background, as dis- 
cussed around (4.7). The scalar profile takes the form 



0(r) = v{r)e 



a{e) 



(B2) 



Now this string carries magnetic fiux F^^ in the bulk, 
but if we are interested in looking at a normalizable state 
then this magnetic fiux must not penetrate all the way 
to the boundary - if it did, it would then define a field- 
theory state with an applied magnetic field, which is not 
what we are looking for. Thus something must happen 
to the magnetic field. It is clear that it must spread out 
along the field theory directions, and thus turn from a 
radial magnetic field into a field parallel to the bound- 
ary, i.e. it turns into a component of F^^ rather than 
F^^. See Figure 8. However this component of the bulk 
field strength tensor can now be interpreted as a circular 
current fiow in the boundary, i.e. 

1 



if 



IQFT 



(B3) 



To demonstrate this, let us compute this current far from 
the vortex core. The dependence on p can be neglected, 
and the relevant Maxwell equation following from (Bl) 
is 

iV^F^') + 2v'qg''V^ {dea - qAe) = (B4) 

This is essentially a forced equation for Ag^ where the 
forcing term is given by 9^ a = 1. We may easily solve it 
by setting. 



Ae{r) = -{l-ae{r)) 



(B5) 



where ae{r) obeys the homogenous part of the same equa- 
tion (B4), i.e. 



-2^r (a 

9f 



-gg^^g^'drae) 



2v\g' 



-gqae = . (B6) 



Furthermore, we require that ao{r) be regular at the hori- 
zon and satisfy the boundary condition ae{r oo) = 1, 
such that at infinity Aq (the field theory source) van- 
ishes. From (B3) is simply given by — ^ times the 
radial derivative of ao at infinity. 

Now note that the equation obeyed by ao{r) is simply 
that for a spatially homogenous gauge fiuctuation on the 
superfiuid background; the presence of the vortex cannot 
be seen in (B6). In fact, since ae is regular at the horizon, 
(B6) is precisely the equation that we solve if we want 
to compute the two point function of the current Gee at 
zero frequency and momentum. It is shown in [52] that 
on such a field configuration the boundary theory Green's 
function is precisely the ratio of the radial derivative drae 
to ae at infinity; converting from Gee to the cartesian Gxx 
we conclude that the circular current satisfies 



je{p) = 



G^^(a; = 0,A: = 0) 



(B7) 



Now in a relativistic superfiuid one way to define the 
superfiuid density ps is via the zero-frequency value of 
e.g. [42-44]: 



G^^(cj = 0,/c = 0) 



(B8) 



We thus conclude that the circulation around the vortex 
is 



dOje 



2nps^ 
q p 



(B9) 



as is expected. 

Now in the calculation in the main text we assume 
that all of the magnetic fiux remains locked in the string 
all the way up to infinity, and thus we are studying the 
field theory with some magnetic fiux inserted into the 
system. We should imagine this as being a pinning po- 
tential for the vortex: the location of this fiux defines the 
position of the vortex ocq in the field theory, and presum- 
ably the force calculated in (4.32) actually acts on the 
apparatus that is producing this fiux. In a truly normal- 
izable vortex state this fiux would not be there, the field 
lines would spread out parallel to the boundary, and the 
vortex would be free to move around. However in this 
situation it is somewhat difficult to precisely define the 
force on the vortex; to extract it we would need to exam- 
ine the equation of motion of a dynamical vortex, which 
appears to be a somewhat more difficult problem then 
the one solved in the text. If such a calculation could be 
done we expect the answer to be roughly the same, as the 
only distinction between these cases is near the boundary 
where the scalar profile (and thus the associated charge 
density) is small. 
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FIG. 8. Two possible vortex solutions. Left: the vortex has magnetic flux enclosed in the string all the way to the boundary. 
This is the vortex that we study in the bulk of the paper, but note that it is not a normalizable solution as it involves a 
nontrivial Fxy source in the field theory. The location of this flux defines the position of the vortex, and the vortex is not free 
to move. Right: the vortex solution analyzed in this appendix. This has no magnetic field perpendicular to the boundary; the 
field tangential to the boundary can be interpreted as circulation, as in (B9). This vortex is a truly normalizable solution, and 
is free to move. 
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